Chiral asymmetry of the spin-wave spectra in ultrathin magnetic films 
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We raise the possibility that the chiral degeneracy of the magnons in ultrathin films can be lifted 
due to the presence of Dzyaloshinskii-Moriya interactions. By using simple symmetry arguments, we 
discuss under which conditions such a chiral asymmetry occurs. We then perform relativistic first 
principles calculations for an Fe monolayer on W(110) and explicitly reveal the asymmetry of the 
spin-wave spectrum in case of wave- vectors parallel to the (001) direction. Furthermore, we quanti- 
tatively interpret our results in terms of a simplified spin-model by using calculated Dzyaloshinskii- 
Moriya vectors. Our theoretical prediction should inspire experiments to explore the asymmetry 
of spin-waves, with a particular emphasis on the possibility to measure the Dzyaloshinskii-Moriya 
interactions in ultrathin films. 
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It is by now well-established that relativistic effects 
play a fundamental role in the magnetism of nanostruc- 
tures, in particular, for thin films and finite deposited 
nanoparticles. Over the past two decades, a vast num- 
ber of experimental and theoretical studies has been 
published to explore related phenomena such as mag- 
netic anisotropics, spin-reorientation phase transitions, 
and non-collinear magnetic orderings. [l|, y, |3|, l4 l5[ 

The antisymmetric exchange interaction between two 



magnetic atoms, Eum = Dij (M; x Mj), where 
and Mj denote the spin-moments of the atoms la- 
beled by i and j, has been proposed 50 years ago by 
Dzyaloshinskii [6(] and Moriya [7j. The Dy is called 
the Dzyaloshinskii-Moriya vector being identical to zero 
if the sites i and j experience inversion symmetry. It 
has been put forward just about ten years ago that 
an enhanced Dzyaloshinskii-Moriya interaction (DMI) at 
surfaces or interfaces can give rise to novel phenom- 
ena in nanomagnetism such as to noncollinear inter- 
layer coupling, [a, [9|1 to unidirectional competing magnetic 
anisotropics, [lOj or to stabilization of non-collinear (chi- 



ral) magnetic orderings. [Ill ll 



A breakthrough on this field happened when the reso- 
lution of spin-polarized scanning tunneling microscopy 
enabled to detect magnetic pattern formation on the 
atomic scale in monolayer-thin films. Such periodic mod- 
ulations have been observed for Mn monolayers deposited 
on W(110) and W(001) and, could successfully be inter- 
preted in terms of a combination of relativistic first prin- 
ciples calculations and a simple micromagnetic model as 
the consequence of large DM interactions. [l3, JJ|- Using 
the same theoretical basis it was even possible to explain 
the homochirality of the domain walls in two monolayers 
of Fe on W(110), [151 in agreement with previous exper- 
imental observation, lfj 

In this Letter, we investigate a consequence of the DM 
interactions on the spin-wave spectra in ultrathin films, 
not yet explored in the literature. We argue that the 
chiral degeneracy of the spin-wave (SW) spectrum can 
be lifted due to the Dzyaloshinskii-Moriya interactions 



and discuss under which conditions such a chiral asym- 
metry occurs. Based on relativistic first principles calcu- 
lations, we explicitly evidence the asymmetry of the SW 
spectrum of an Fe monolayer on W(110) in case of wave- 
vectors parallel to the (001) axis. We then quantitatively 
interpret our results in terms of a simplified spin-model 
by using calculated DM vectors. By emphasizing the pos- 
sibility of probing the DM interactions in ultrathin films, 
we impel experiments to explore the proposed effect. 

We start our study with simple considerations based 
on classical spin-waves. If the atomic magnetic mo- 
ment in the ground state of a ferromagnetic monolayer is 
Mo = Moeo with eo being a unit vector, then a spin- wave 
of wave-vector q and chirality index (rotational sense) 
c = ±1 is defined by the magnetic orientations, ei(q, c) = 
ni cos(qRi) sin0 + cri2 sin(qRi) sin0 + eocos0, where 
ni _L eo and ri2 = ni x eo are unit vectors, is the 
position vector of site i and the is the relative angle 
between the moments and eo. Inspecting the energy of 
the SW in terms of an extended Heisenberg model con- 
taining tensorial exchange interactions. (l7j it turns out 
that only the antisymmetric exchange interactions, i.e., 
the DM interactions give rise to a chirality dependent 
contribution, 

E D m(cl, c) = csin 2 0^ (D y - • e ) sin (q(Rj - R,-)) . 

(1) 

The above expression also implies that only the compo- 
nents of the DM vectors parallel to eo influence the SW 
energy and that a reversed chirality can be converted 
into a propagation of the SW in the opposite direction: 
Edm(<1, -c) = E DM (-q,c) = -E DM (q,c). 

The orientations of the DM vectors in a ferromagnetic 
monolayer have been analyzed for different 2D lattices 
in Refs. 11 1 and [lij]. In particular, if the lattice has a 



twofold rotational symmetry, such as in case of the (001) 
and (011) surfaces of cubic lattices, all the DM vectors 
lie in-plane. Clearly from Eq. fl}, a chiral asymmetry 
of the SW occurs then only for an in-plane ground-state 



2 



magnetization. Furthermore, if the ground-state magne- 
tization is in a mirror plane of the monolayer no chiral 
asymmetry applies for wave-vectors being parallel with 
e . 

In order to demonstrate the chiral asymmetry of the 
SW's we have chosen a ferromagnetic Fe monolayer de- 
posited on W(110), since (i) it exhibits an in-plane 
ground-state magnetization 1!J 2(| and (ii), as for the Fe 
double-layer flEI or for a Mn monolayer on W(110) and 
W(001), [H,H[Il| large DM interactions are expected. 
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FIG. 1: (Color online) Sketch of the lattice positions (left) 
and of the surface Brillouin zone (right) of a bcc(llO) plane. 
The high symmetry points of the surface Brillouin zone are 
also labeled. 

Our notation used for the principle axes of a bcc(llO) 
plane are shown in Fig. [1] depicting the structure of the 
real-lattice and the surface Brillouin zone with the high 
symmetry points. The theoretical bec W lattice con- 
stant as given in Ref. [22|, a = 3.205 A, was chosen for the 
in-plane lattice constant (along the Y axis) throughout 
the system. All the interlayer distances were fixed to the 
ideal bcc(110) value, d = s/2a/2 = 2.266 A, but the inter- 
layer distance between the Fe and the topmost W layer 
was relaxed by -12.9 % (d Fe - W = 1.974 A) according 
both to experiment [23| and to theory [24|. The calcu- 
lations were performed in terms of the fully relativistic 
Screened Korringa-Kohn-Rostoker (SKKR) method [25j 
by using the local density approximation and the atomic 
sphere approximation (ASA). It should be noted that the 
SKKR method makes use of a semi-infinite geometry for 
the substrate, therefore, the calculations are not affected 
by ambiguities related to a supercell or film geometry. 

We calculated magnetic anisotropy energies, E(00l) — 
E(lT0) = 2.11 mcV and £(110) - £(110) = 0.41 mcV. 
This implies that, in agreement with other theoretical 
works 126 



27 1 and with the experiment, [l9(, the ground- 
state magnetization of FeW(llO) is in-plane with an easy 
axis along the (110) direction and the hard axis is along 
the (001) direction. It should be noted that the magne- 
tostatic dipole-dipole interaction also favors the (110) di- 
rection by about 0.01 meV with respect to the (001) axis 
and by 0.11 mcV with respect to the (110) direction. [28j| 
We applied a recent relativistic extension [17| of the 
torque method [2§| to evaluate tensorial exchange inter- 
actions for FeW(llO) from first principles. This method 
opened the way to atomistic spin-model simulations of 



nanostructures accounting for relevant relativistic inter- 
actions, such as the on-site magnetic anisotropy, the 
anisotropic symmetric exchange interaction and the an- 
tisymmetric exchange interaction. [Til. [21 . 

By using the convention, H = — Jij e i e ji our cal- 
culated isotropic exchange interactions for the first few 
neighbors are, J i = 10.84 meV, J02 = —3.34 meV, 
J03 = 3.64 meV and J04 = 4.60 meV. Note that, in par- 
ticular, the nearest neighbor interaction, J01, is about 
four times less in magnitude than the corresponding pa- 
rameter in Ref. 27. To lend confidence to our values 



for Jy, we performed Monte Carlo simulations and ob- 
tained a Curie temperature of about 270 K, in very good 
agreement with experiment (225 K).[l|| Note that ran- 
dom phase approximation (RPA) calculations in Ref. [13 
provided a Tq above 1000 K which is, most likely, the 
consequence of the overestimated NN exchange interac- 
tion. 

By using a canonical quantization of the linearized 
Landau-Lifshitz equations, we also developed a method 
to calculate the adiabatic SW spectra of bulk and lay- 
ered systems on a relativistic first principles basis. 17 
Although, within this approach, the interaction of the 
spin-waves with the Stoner continuum is neglected, the 
main features of the SW spectra due to relativistic effects 
are expected to be well-described. Notably, in case of a 
monolayer, two SW solutions are obtained with the ener- 
gies, E + (q) and E~(q), that correspond to the chirality 
indices +1 and -1, respectively. 
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FIG. 2: (Color online) Calculated adiabatic spin-wave spec- 
trum of FeW(llO) along the X axis as given in Fig. [TJ 

In Fig. [2] the calculated adiabatic SW spectrum is 
shown along the X axis. It was demonstrated in Ref. [27 
that the adiabatic SW energies and the SW dispersion 
obtained from RPA agree well for wave-numbers as large 
as about 1 A -1 . We, therefore, display the adiabatic SW 
spectrum for only \q\ < 1.2 A -1 . Anticipated from the 
symmetry analysis above, since in this case the ground- 
state magnetization and the wave-vectors lie in a mir- 
ror plane of the system, the spectrum is degenerate, i.e., 
E + {q) = E~(q). Correspondingly, the SW dispersion is 
symmetric, E(q) = E(—q). Note that the energy range 



3 



of the disp ersion in Fig. [5] is approximately half of that 
in Rcf. |27|, which we again attribute to the very different 
exchange interaction parameters in the two theoretical 
works. 
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FIG. 3: (Color online) Calculated adiabatic spin-wave spec- 
tra with chirality index +1 (triangles) and -1 (spheres) of 
FeW(llO) along the Y axis, see Fig. Q] 



Next we inspect the SW spectrum for wave-vectors 
parallel to the (001) axis displayed in Fig. [3] Since in 
this case q is perpendicular to the ground-state mag- 
netization, our symmetry analysis predicts lifting of the 
chiral degeneracy of the spectrum, which can evidently 
be inferred from Fig. [3] Furthermore, the relationship 
E + (—q) = E~{q) is clearly regained. As compared with 
Ref . [27|, again a difference by a factor of two in the energy 
range of the magnons can be noticed. 




FIG. 4: (Color online) Sketch of the calculated 
Dzyaloshinskii-Moriya vectors between an atom (C) and its 
nearest (1) and next nearest (2) neighbors in an Fe monolayer 
on W(110). 

In order to demonstrate that the observed asymme- 
try of the SW spectrum results from the DM interac- 
tions we performed a model calculation for AE(q) = 
E + (q) — E~(q). Our first principles calculations indi- 
cated that the DM vectors for the nearest and next near- 
est neighbors, visualized in Fig. 01 arc at least by an order 
larger in magnitude than the ones for more distant pairs. 
By using the method described in Rcf. [17| the asymmetry 



of the SW energy can then be expressed as, 

AE(q) = ^Df sxn(ha) - ^£>f sin(ga) , (2) 



A/n 



Mn 



where M = 2.22 [i^ is the spin-magnetic moment per 
atom and Df and Df are the magnitudes of the x com- 
ponents (parallel to the ground-state magnetization) of 
the DM vectors for the nearest and second nearest neigh- 
bors, respectively. 
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FIG. 5: (Color online) Squares: asymmetry of the spin-wave 
spectrum of FeW(llO) along the Y axis as derived from the 
data in Fig. [3j solid line: the function, Eq. ([2}, obtained from 
a second nearest neighbor model with the calculated DM in- 
teractions. 

In Fig. [5] we plotted the asymmetry, AE(q) 1 of the SW 
spectrum of FeW(llO) along the Y axis obtained from 
the data in Fig. [3] As can be inferred from this figure, 
AE(q) exhibits local extrema at about q = ±0.44 A -1 
with \AE(q)\ ~ 15 meV and changes sign at q = ±0.83 
A -1 . Apparently, these features of AE(q) are fairly well 
reproduced by the function, Eq. ([2]), when using the cal- 
culated parameters, Df = 1.42 meV and £>f = 6.08 meV. 
Note that the characteristic extrema of AE(q) are deter- 
mined by the DM interactions between the next nearest 
neighbors, since the sin(ga) function in the second term 
on the rhs of Eq. ([2]) reaches a maximum/minimum at 
\q\ = it /2a ~ 0.49 A -1 . The deviations of the asymmetry 
of the SW energy from this model function are related 
to the DM interactions between more distant pairs that 
add low-frequency modulations to the SW dispersion. 

The magnon spectrum of FcW(llO) along the Y 
direction has been measured very recently by us- 
ing spin-polarized electron energy loss spectroscopy 
(SPEELS),[32j a highly suitable technique to probe high 
wave-vector magnetic excitations of ultrathin films. Sur- 
prisingly, the measured magnon energies are about half of 
the theoretical values reported here and smaller by even 
a factor of four than the calculated values in Ref. |2jJ We 
are, however, aware of linear response calculations [33| 
that provided with a very similar magnon dispersion 
along the Y axis as compared to that in Fig. [3] Thus, the 
low energy of the measured magnon spectrum [32| should 
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most probably be attributed to effects not included in 
the fi rst p rinciples calculations, such as s pin -charge cou- 
pling [34| or phonon-magnon interaction. [35[ 

Considering the size of the SW asymmetry obtained 
from our calculations, e.g., about 20 % at q = ±0.44 
A -1 with respect to the average energy, (E + (q) + 
E~(q))/2, we strongly suggest that it should be acces- 
sible to experiments. Indeed, preliminary measurements 
on FeW(llO) [36[ indicate the presence of an asymmetry 
in the magnon spectrum being quite similar in size and 
shape as in Fig. [5J Further candidates for experimental 
observation of the proposed SW asymmetry are ferro- 
magnetic monolayers on substrates with large spin-orbit 
coupling and polarizability (W, Pt or Ir). In case of an 
out-of-plane ground-state magnetization, a small mag- 
netic field might be applied to orient the magnetization 
in-plane, in order to fulfill the necessary condition for the 
chiral asymmetry of magnons. Our model calculation, see 
Eq. @ , clearly implies that such experiments would serve 
as a unique tool to measure the Dzyaloshinskii-Moriya in- 
teractions in ultrathin ferromagnetic films, to be directly 
compared with the results of ab initio calculations. Con- 
cerning, in particular, the role of relativistic effects, such 
a progress would clearly assist a deeper understanding of 
the magnetism in nanostructures. 
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